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Abstract 


Neutrino induced deep inelastic charm production off strange sea partons in the nucleon 
is studied in the framework of the MS fixed flavor factorization scheme. The momentum 
(z) distributions of the produced D mesons as calculated in the next-to-leading order 
QCD analysis are presented and compared to their leading order counterparts. Pertur¬ 
bative stability within this formalism is demonstrated and the compatibility of recent 
next-to-leading order (NLO) strange quark distributions with available dimuon data is 
investigated. Our NLO (MS) results provide, for the first time, the tools required for a 
consistent and complete NLO (MS) analysis of charged current dimuon data in order to 
extract the strange sea density in NLO. It is furthermore shown that the absolute predic¬ 
tions of radiatively (dynamically) generated strange sea densities are in agreement with 
all present measurements. 



Recent attempts 00 to determine the strange sea content of the nucleon , s(x, Q 2 ), 
from a next-to-leading order (NLO) QCD analysis of neutrino induced charm production 
(opposite-sign dimuon events) were undertaken in the framework of the NLO formalism of 
0| . In this formalism one considers, in addition to the leading order (LO) process W + s —> 
c, etc., the W + g —> cs gluon fusion contribution with m s 7 ^ 0, i.e. one employs a finite 
mass regularization and subtracts from it that (collincar) part which is already contained 
in the renormalized, Q 2 -dependent s(x, Q 2 ). In a previous publication 0 we have started 
a study of NLO QCD effects based on the conventional MS factorization scheme where the 
mass of the strange quark is neglected ( m s = 0 ) and the corresponding mass singularities 
are dimensionally regularized [h]. Such an analysis is imperative in connection with all 
those NLO parton distributions and cross sections given in this, most frequently utilized, 
massless MS factorization scheme. The analysis in 0 was incomplete, however, in the 
sense of not addressing the detailed momentum distributions of the produced c-quark 
(D-meson) utilized in 0, 0. The present paper is devoted to this issue and the relevant 
expressions are now calculated in their full differential form, by extending the original 
inclusive analysis in 0 . 

The formulae and notations of 0 will be strictly taken over everywhere, being extended 
merely by the introduction of the fractional charm momentum variable ( = p c • Pn/q ■ Pn, 
i.e. H 9 ' 9 , p 2 , X) in 0 will be replaced by H?' 9 (£',(, p 2 , A) where the fermionic and 
gluonic NLO coefficient functions Hf =1 2 3 and H 9 for heavy quark (charm) production 
derive from the subprocess W + s —> eg (together with the virtual vertex correction to 
W + s —> c) and W + g —> cs, respectively. Furthermore, d 2 a u X) / dxdy in 0 will be re¬ 
placed by d 3 cr u ^/dxdydz with the relevant semi-inclusive structure functions, describing 
the momentum distributions of the produced D-mesons, being obtained from a convolu¬ 
tion with the charm fragmentation function D c (z), [] 

Rnstead of z = po ■ Pn/q ■ Pn, it is customary to use |[j, [t| z p = pd/Pd^ which differs only 

marginally from z in the relevant experimental region (z 0.05). 
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where Ff = F, c , Ff = F 3 c /2, F 2 = F 2 c /2£ with £ = x(l + m 2 c /Q 2 ), Cmm = "hr/s = 
(1 — A)£'/(l — A£') and A = Q 2 / [Q 2 + ml). The H 9 ’ 9 are given in the Appendix and our 
choice for the factorization scale will be fi 2 F = Q 2 + ml, although the results are not very 
sensitive to this specific choice. The fractional momentum of the D-meson is denoted by 
z = p D ■Pn/Q'Pn and the (factorization scale independent) charm fragmentation function 
is taken as 0 

D c (z) = N jz 1 - z~ 1 - e c /(l - z) | (2) 

with the normalization constant N being related to e c via /q 1 dzD c (z ) = 1, while e c 
itself will be fitted to the measured z-distributions. One can infer these distributions by 
utilizing the LO version of eq. (1) which has a simple factorized form 


d 3 a u L g> dvr (p) 


' LO 


D c (z) 


dx dy dz dx dy ’ 
together with the corresponding parameters in the LO analysis 0 of which = 0.20 ± 
0.04 is the most relevant now @j. It should be mentioned that the CCFR group has 
also performed a ’NLO’ analysis JI[] of their dimuon data using the simple factorized LO 
expression in eq. (3) instead of the correct (convoluted) NLO expressions in (1), which are 
just the finite-mass extensions of the well known NLO results for light (massless) quarks 
0. We shall analyze the ’data’ [0 as parametrized in LO by eq. (3). 


To simplify the presentation of our results we define Q 

- 1 < l + Q 2 ! M w) 2 ir 
£s(t,z,Q) eff - 2 G 2 fMnEv I v * 


1-2 




(4) 


dx dy dz 

which has also been studied experimentally |lj and where the superscipt cs refers just to 
the CKM non-suppressed (VA) component of s' = |VA|"s + \Vcd\~ {d + u)/2 in eq. ( 1 ). 
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In LO the cross section in (4) reduces to [||, cf. eqs. (1) and (3), 


£s(£,z,Q 2 ) eff = (1 


m: 


^(£,/i 2 ) D c (z) 


( 5 ) 


2 M^E U C, 

Our LO and NLO results for s e ff are shown in figs. 1 and 2 using the GRV [EilJ and 
CTEQ4 [ilTJ parton densities, respectively. It should be emphasized that the ’data’ on 
d 3 a/dxdydz are well reproduced in NLO provided one selects a harder nonperturbative 
charm fragmentation function D c (z ) corresponding to £^ LO = 0.06 in eq. (2) as expected 
to be necessary for compensating the softening effects due to the NLO 0(a s ) corrections 
(as is evident from comparing the dashed and solid curves in figs. 1 and 2). Similar 0(a s ) 

-NLO _ 


softening effects due to gluon bremsstrahlung in e + e 


DX imply 
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0.06 ± 0.03. It is interesting to note that similar values for e ^ LO , i.e. £^ LO — 0.06 are 
obtained in NLO by comparing the predictions of the NLO perturbative fixed order QCD 
corrections to e + e _ —> DX, based on the corresponding cross sections in ]Li, with the 
(v/i~ 10 GeV). 


ARGUS data 


Apart from the very large z region (z ~ 0.8), where the perturbative ln(l — z) singu¬ 
larities have to be resummed [|1^] and nonperturbative higher twist contributions ]17| will 
play a dominant role, our LO and NLO results in figs. 1 and 2 show a remarkable pertur¬ 
bative stability throughout the experimentally relevant z and x region (z ~ 0 . 8 , 0.01 ~ 
x ~ 0.4) for both sets of parton densities, despite the fact that we have (arbitrarily) fixed 
m c = 1.5 GeV in LO and NLO. It should be reemphasized that a proper utilization of the 
corresponding LO and NLO parton densities and coupling constants (and possibly m c ) is 
essential for the observed perturbative stability of the physical cross sections in eqs. ( 1 ), 
(4) and (5). The results in Eg. 1 support the dynamically predicted s(x,Q 2 ), generated 
purely radiatively from s(x,yU 2 ) = 0 at ti 2 LO nlo = 0.23, 0.34 GeV 2 [[L0j. On the other 
hand, the CTEQ4 densities overshoot somewhat the shaded ’data’ band due to the 
larger strange sea densities. 


The observed perturbative stability, as shown for example in fig. 1, can possibly be 
improved by adopting a different charm mass in LO and NLO, e.g. m c = 1.3 GeV and 
m c = 1.7 GeV, respectively, following refs. 10 - The observed variations of the LO and 
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total NLO results are shown in fig. 3. 


In order to determine the NLO strange sea density, the CCFR collaboration a, a 
employed, as discussed at the beginning, the NLO mass-regulated (m s ^ 0) W + g —> cs 
gluon fusion contribution within the ACOT framework || where the NLO 0(a s ) quark- 
initiated contribution (W + s —► eg) is disregarded (since it has not been calculated yet 
for m s 7 ^ 0). An explicit (more differential) calculation of this gluon-initiated production 
process along the lines of H shows [I 8 [] that cl 3 cr/dxdydz becomes insensitive to rn s pro¬ 
vided one chooses rn s ~ 200 MeV, as is usually done 00 and practically coincides with 
our MS results, e.g. the dashed-dotted curves in fig. 1. Moreover, these results depend 
very little on the specific choice of the initial gluon density for the experimentally relevant 
(x, Q 2 ) domain, as is obvious from the dashed-dotted curves in figs. 1 and 2 ; the gluon 
density, obtained in the original CCFR analysis 00- gives again a very similar result. It 
should, however, be noted that CCFR fl] approximated the full convolutions of the NLO 
gluonic contribution in eq. (1) by a factorized expression similar to the one in eq. (3) 
while keeping the LO D c (z ) with e^° ~ 0.20. This ’approximated’ factorized expression 
yields rather different and ^-distributions. Despite the fact that such a factorized NLO 
approach is theoretically not justified, this inconsistency could be partly ’compensated’ 
by choosing a different value for e c (and possibly also for m c ). In the latter case the 
extraction of s^lo becomes rather ambiguous, since the utilized value of e c cannot be 
consistently obtained within an analogous NLO analysis of, say, e + e _ —> DX. On top of 
this comes the contribution from the NLO quark-initiated subprocess W + s —> eg and the 
NLO vertex correction to W + s —> c, so far neglected by the CCFR analysis [jf], [2J, which 
is not negligible according to the dotted curves in figs. 1 and 2. It is thus clear that a 
consistent and complete NLO-reanalysis of the CCFR dimuon data would certainly be 
worthwhile. 

For illustration we finally compare the NLO strange sea density inferred by the CCFR 
collaboration [jTJ with their LO one in fig. 4, where the GRV94 and CTEQ4 strange den¬ 
sities are also shown at Q 2 = 4 and 22 GeV 2 . The three different strange sea densities are 
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comparable in the limited x-domain of the CCFR dimuon data at fixed Q 2 (as indicated). 


Outside this region the enhancement of the ’extracted’ snlo{x,Q 2 ) with respect to, say, 
the GRV sjstlo should not be taken too literally since the extraction beyond the measured 
region depends on the assumed input x-shape as well as on the Q 2 -evolution.| On the 
other hand the excess of the LO/NLO CTEQ4 predictions over the CCFR ’data’ band 
(as extracted in LO) in Eg. 2 is caused by the strong enhancement of the CTEQ4 strange 
quark densities with respect to the LO/NLO GRV ones, in particular for x ~ 0.03. 

To summarize, we have calculated and analyzed within the MS fixed flavor scheme the 
momentum (z) distributions d 3 a/dxdydz of the D mesons produced in neutrino induced 
deep inelastic charm production off strange sea partons (opposite-sign dimuon events) 
according to the quark- and gluon-initiated NLO subprocesses W + s —> eg (together 
with the vertex correction to the LO Born process W + s —> c) and W + g —> cs, respec¬ 
tively. Perturbative LO/NLO stability within this formalism has been obtained in the safe 
0 ~ 0.8 region and the compatibility of recent LO and NLO strange quark sea densities 
with available CCFR dimuon data §0 has been investigated. We studied also the mass- 
regulated (m s 7 ^ 0) W + g —> cs gluon fusion contribution within the ACOT framework || , 
originally used by CCFR for extracting a NLO strange sea density Jl]| from their dimuon 
data, where the NLO 0(a s ) quark-initiated W + s —> eg contribution is disregarded. This 
m s 0 gluon initiated contribution turns out to be insensitive to m s provided one chooses, 
as usual, m s ~ 200 MeV and practically coincides with our MS results [0 . Nevertheless 
the CCFR NLO analysis |[[ of the dimuon data for extracting snlo(x, Q 2 ) is inconsistent 
due to the use of a theoretically unjustified ’approximated’ factorized expression in x 
and z for d?o v ^ v ’/dxdydz, as well as incomplete due to the neglect of the quark-initiated 
0(a s ) contribution. A consistent and complete NLO-reanalysis of the CCFR dimuon 

2 It should be noted that CCFR has assumed a constant value for 


k(Q 2 ) = [xs(x, Q 2 ) + xs{x, Q 2 )] dx \xu(x, Q 2 ) + xd(x, Q 2 )] dx 




using k(Q 2 ) ~ 0.48. This is partly responsible for an enhancement of s(x,Q 2 ) at smaller values of Q 2 , 
since a LO/NLO QCD Q 2 -evolution gives k (4 GeV 2 )/ft(20 GeV 2 ) ~ 0.8. 
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data, using the fully differential and complete NLO (MS) contributions as presented in 
this article, would therefore be worthwhile. 
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Appendix 


The fermionic NLO (MS) coefficient functions Hf for heavy quark (charm) production in 
eq. (1), calculated from the subprocess W + s —»■ gc and the vertex correction to W + s —> c 
are given by 


HUSX.AA) = 3(1-C) ( 

l l l F 


1 - £ + (1 - 0 In £ 91 - 


£(i - A 0 


1 -^ 


1 in.^ 1 -^ 


i-e i-Ae 


+ - |-5(1-05(1-0 


1 / 1 + 3A ^ l\ ^ 7T 2 

K a + T + 4 + — 
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1 e +(1-0^ Ae 


( 1-0 


2 e 


(i- 0 + (i-O 

1 — A£ 


o-e/ 

i-(i-0 


A 

1 -e 

(1 - A 0 2 

1 -A^ 


i-£ 


i-(i-0 


i-e 


(Al) 


#f(0 0/4, A) = ^(OC,^,A) + 3|<5(l-05(l-0^ 

- 2^(1-3A)[1-(1-C)\^] + (1-A) N 

#1(00/4,A) = #?(00/4,A) H- 2 |/(1 - 0[1 — (1 - C) 


1-Ag 

1 -^ 


] - (1 - A 0 


(A2) 

(A3) 


1 4 (\ i t2\ 

with = —(1 —A) ln(l —A) , P®(£) = - -— and where the distributions 

A ' 3 V 1 - £ 7 4 - 

are defined by 


^ /( 0 0 ^/( 0 -/(I) 

3 ^ (1-0+“0 ^ l-e " 


r * a 

' C min V J- S> 


/' + /(c) _ ( (1) < A4 > 
' Cmin -i- S> 


with <C m j n = (1 — A)£/(l — A£). When integrated over £, these results reduce to the final 
inclusive results obtained in [[|, i.e. 


r <K + , «,c, 4 +) 

** Lmin 


(AS) 
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The gluonic NLO (MS) coefficient functions Hf for heavy quark (charm) production in 
eq. (1), as calculated from the subprocess W + g —> cs, are given by 


h‘((, 


99 

1 1 

c 


In 


Q 2 + m 2 c 

0 F 


+ In 


(i-0 s 


( 1-0 


e(l - AO 

+ fcftec.A) 


+ «!-«) 


(A 6 ) 


where 


^?(0C,A) = -^(l-A)(l-2A) + ^(l-A)(l-2AO 

+ 2fA(l - AO - 1 (A7) 

^f(0C) A) = ^(1 — A)(l — 6A + 6A 2 ) H—^-(1 — A)(l — 2A0 

+ A[6A£(1 - AO - 1] (AS) 

m, c, A) = 2(1 - A)(l - 2 A) - |[PW(0 + 01 - e - A0(1 - A)] (A9) 

with 0^(0 = -[0 + (1 — 0 2 ] and the © distribution is defined in (A4). When 
integrated over 0 these results reduce to the final inclusive results obtained in J|, |J, i.e. 

mt,i4,\) = [ 1 dcm 00 a) 


(A10) 
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Figure Captions 


Fig. 1 The predicted £s(£, z,Q 2 ) e ff in eqs. (4) and (5) as calculated in LO [NLO] using 
the dynamical (GRV) parton distributions of ref. [[TIJ and the charm fragmentation 
functions in eq. (2) corresponding to e c = 0.20 [0.06]. The individual NLO quark- 
and gluon-initiated contributions derive in an obvious way from eq. (1). The charm 
mass has been (arbitrarily) fixed in LO and NLO, m c = 1.5 GeV. The shaded band 
represents the ’data’ as extracted from [0] via eq. (3). 


Fig. 2 Same as fig. 1 but using the LO CTEQ4L and the NLO (MS) CTEQ4M parton 
densities 0 - 


Fig. 3 LO and total NLO predictions using different values for m c (thin curves). The LO 
and NLO results for m c = 1.5 GeV (thick curves) are the same as the corresponding 
ones in fig. 1. The predictions are based on the GRV densities [10]. The shaded 
experimental ’data’ band corresponds, however, to m c = 1.3 GeV as utilized in LO 
in 0,0 and is the same as in fig. 1 . 


Fig. 4 Strange sea densities in LO and NLO at Q 2 = 4 GeV 2 and Q 2 = 22 GeV 2 
according to CCFR |]J and CTEQ4 |lT|. The GRV strange densities |L(J refer to 
absolute QCD predictions since they have been generated purely radiatively from 
a vanishing input s(x, /i 2 LO NLO ) = 0. The NLO CCFR strange density corresponds 
to the choice of the factorization scale Q 2 + m 2 in Table 5 of ref. [T]. The actual 
x-ranges covered by the CCFR data [7] at these fixed values of Q 2 are indicated. 
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